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The high-energy factorization and the associated B-JIMWLK or BK evolution equations 
are presented, using the example of DIS structure functions. The necessity of taking gluon 
saturation into account is discussed, and also the various approximations underlying high- 
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1. Introduction 

QCD hard processes are described within the improved parton model, thanks to 

the collinear factorization. In that kinematical regime, only a few quasi on-shell 

quarks and gluons per hadron are relevant for the dynamics, and thus hadrons are 

appropriately described by parton distribution functions. 

However, this picture breaks down in the high-energy limit or equivalently at 

low Bjorken-x. Indeed, due to the enhancement of soft gluon Bremsstrahlung, the 

gluon occupation number reaches non-perturbatively large values at low enough x. 

Gluons in this regime are better described as a semiclassical field rather than as a 

dilute gas of partons. The growth of the gluon occupation number is then tamed 

by the no n-linearity of the Yang-Mills equations. This is the phenomenon of gluon 
saturation.HEEEE] 

This phenomenon is taken into account in various formalisms which aim at 
describing dense-dilute and/or dense-dense observables. In the former case, the low- 
x gluons of a hadron or nucleus target are probed by a dilute projectile. Typical 
examples of that case are Deep Inelastic Scattering (DIS) observables at low- a; and 
forward physics in hadron-hadron or hadron-nucleus collisions. By contrast, dense- 
dense observables involve the collision of the low-cc gluon field of two hadrons or 
nuclei against each other. The inclusive n-particle production at mid-rapidity in 
heavy ion or hadron-hadron collision at high-energy enters in that category, as well 
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as the n-point correlators of the energy-momentum released at mid-rapidity in such 
collisions. 

In this contribution, I focus on the paradigmatic example of DIS, whereas 
dense-dense processes are discussed by Venugopalan in his contribution.® Note 
that both for dense-dilute and dense-dense observa bles, the high-energy l eading 
logs are resummed by the same B- JIMWLK equation s] ' ' ' ' ' ' ' ' ' acting 
on Wilson line operators, first derived without explicit reference to the projec- 
tile. Hence, there is a factorization property valid at leading logarithmic accuracy 
for several very different observables. In the large N c limit, the B- JIMWLK equa- 
tions ap plied to a dipole of Wilson lines operators reduce to the much simpler BK 

equation! 7 | 1 6 1 1 7| B oth the B - JIMWLK equations and the BK equation reduce to 
11 811 QI90I 

the BFKL equatio n ° | a | U | in the dilute limit for the t arget. 

Recently, Balitsky and Chirilli have calculate d 2 ^ 22 ! the evolution equation re- 
summing the next to leading logarithms (NLL) in the particular case of dipoles of 
Wilson lines operators (and thus the NLL BK equation) , as well as the NLO impact 
factor 2 -^ (or high-energy coefficient-functions) for DIS structure functions. These 
calculations show that the high-energy factorization remains valid at NLO/NLL 
accuracy, at least in the case of DIS stru ctu re functions. This high-energy factor- 
ization generalizes the dipole factorization^!, valid at LO only. Hence, high-energy 
factorization is now a good candidate for a factorization property valid to all orders, 
analog to collinear and TMD factorizations. 

However, the NLL BK equation is plagued by problems inherited from its lin- 
earized analog, the NLL BFKL equationPSESlln both equations, the NLO contribu- 
tions to the kernel are large, especially in the collinear and anti-collinear kinematical 
regimes, signaling a breakdown of the weak coupling expansion in most of the avail- 
able phase-space. Hence, in order to obtain stable and reliable solutions of the NLL 
evolution equatio ns, one ne eds to perform appropriate resummations, already done 
in the BFKL case! 27 ! 28 ! 29 ] 

In the next section, the derivation of the high-energy factorization is outlined 
in the context of light-front quantization, and the underlying assumptions are dis- 
cussed in detail. In the section [3l the large corrections appearing in NLL BFKL 
and NLL BK equations and their physical origin are discussed. Progresses towards 
the resummation of these corrections to the NLL BK equation are also reviewed. In 
the section [4j preliminary results about the resummations of the large corrections 
of kinematical origin are introduced, with emphasis on the factorization scheme 
dependence of these effects. 



2. High-Energy Factorization for DIS 

The high-energy factorization for DIS can be obtained both from covariant or light- 
front approaches to QCD. The formali sm use d by Balitsky and Chirilli and dis- 
cussed in their respective contributions J2012U j s fully covariant and thus suitable 
to perform systematically higher order calculations. As an alternative, let us use a 
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light-front formalism, which allows to make the physical content of the calculations 
more transparent at leading order or at leading logarithmic accuracy. This formal- 
ism is essentially the QCD version of the one constructed by Bjorken, Kogut and 
Soper for QED.El 

Let us choose a frame such that the projectile moves mainly along the x + direc- 
tion and the target mainly along x~ ', and such that the target is much more boosted 
than the projectile. In the high-energy (or low x) limit, the target is more and more 
Lorentz contracted along x + , so that the x + time interval around x + — needed 
for the projectile to pass through the target becomes arbitrarily short. Hence, the 
forward scattering process related to DIS by the optical theorem happens in three 
steps 

• the projectile propagates in the vacuum form x + = — oo to 0, 

• the projectile interacts with the target at x + = 0, 

• the projectile propagates in the vacuum form x + = to +oo, 

up to contributions suppressed by powers of the Lorentz 7 factor of the target. 

The formalism of light-front wave functions allows to describe the state of the 
projectile just before or just after the collision. A light-front wave function is indeed 
the overlap of the physical state of the projectile with a given Fock state at x + =0. 
In the DIS case, the projectile is in principle the incoming lepton, but as in covariant 
calculations, one can factor out the leptonic tensor, and take effectively the projec- 
tile to be the exchanged virtual photon, which can be transverse or longitudinal. 
At LO, only the quark-antiquark Fock component of the virtual photon projectile 
contributes to DIS, so that the hadronic contribution to the projectile state writes 
formally 

|Proj) L ° = J dPS m *lo )qq) ( (1) 

with a sum over the phase space and quantum numbers of the quark and antiquark. 
At NLO, the quark-antiquark-gluon Fock component contributes as well, so that 

|Proj) NL ° = J dPS m V^NLO |q - } + J dpSmg ^NLO |q - g) (2) 

In the high-energy limit, the various partons present in a Fock state of the 
projectile scatter independently on the target, up to power suppressed corrections. 
Indeed, the time interval for the interaction between the target and the projectile 
goes to zero when boosting more and more the target, and also the typical time 
scale for interactions within the projectile is Lorentz dilated when boosting the 
projectile. 

Following the Color Glass Condensate (CGC) form 

a l isn J3 | 4 | 5 | 13 | 14 | 15 | 

let us 

model the target by a random classical gluon field, obeying some classical statistical 
distribution. I will discuss the validity of this assumption later on. Then, we need 
to know how a parton scatters in the high-energy limit on a classical gluon field 
shrinking to a Shockwave at x + — 0. When passing through such a Shockwave, the 
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transverse position of a parton has no time to change. Hence, it is convenient to 
write the Fock state basis in a mixed representation where the transverse position 
rather than the transverse momentum of the partons is specified, in addition to 
their light-cone momentum k + . Both the probability of helicity flip of the parton 
and the probability of loss of a non-negligible fraction of his light-cone momentum 
by scattering on the gluon Shockwave, are power suppressed in the high-energy 
limit. Hence, the scattering of independent partons on the classical gluon Shockwave 
reduces to the so-called eikonal scattering, which is diagonal in the Fock state basis 
except with respect to color indices. A parton is thus only color rotated when passing 
through the Shockwave field. This color rotation is given by a Wilson line in the 
appropriate representation, defined along the trajectory of the parton through the 
target gluon field, with fixed transverse position x, noted for example U(x) in the 
case of a quark. 

In particular, the S-matrix for the forward scattering of a free quark-antiquark 
color singlet dipole (at transverse positions x and xi) on a classical gluon shock- 
wave is the identity matrix in Fock space times the factor 



1 



S X0X1 = — tr (C/(x )[/t( Xl) ) 



(3) 



Using Fierz identities, the analog S-matrix factor for a quark-antiquark-gluon color 
singlet can be written as 



1 



^x xix 2 _ I 



tr (t/(x )[/t( X2 )) tr(f/(x 2 )f/t( Xl )) tr(C/(x )C/ t (x 1 )) 



where x 2 is the transverse position of the gluon, so that 



<Sx XiX2 



£x 2 xi 



AT2-1 

All in all, the projectile-target total cross section writes at LO as 

|2 



at, =2 J dPS m |^0| 



l-(Sx 0Xl ) 



(4) 
(5) 

(6) 



with a different wave function whether the virtual photon projectile is trans- 
verse or longitudinal, and where ( • • • ) represents the statistical average over the 
target gluon field. And at NLO, one expects the total cross section to be 



2 / dPS qq |* qq 



NLO 1 2 



l-(Sx 0Xl ) 



1 (5 , x xix 2 ) 



(7) 



and so on at higher orders. Using the probability conservation relation 



NLO | 2 

qq 



+ / dPS z * 



NLO 1 2 

qqg I 



l* LO l 
I qq I 



(8) 
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the total cross section at NLO can be rewritten as 



2 / dP5* aa * 



2 



LO 

qq | * qq 



2 



1 (S'xoxi) 



2AT. 



NLO 1 2 



iV? — 1 J qqs 1 ^ 



(^X Xi) (5'xoX2 5x2X1 ) ) (9) 



separating the LO and NLO contributions from each other. 

However, this is not the end of the story. Indeed, the operators constructed from 
light-like Wilson lines such as 5 Xo xi suffer from rapidity divergences, coming from 
gluons of large k + and/or vanishing k~ belonging to the field of target. And the 
usual logarithmic divergence of soft gluon Bremsstrahlung shows up in the NLO 
contribution as 



|2 



dk+ 



dPSqqg |*qq g u | oc — for small k + of the gluon. (10) 

These two divergences are both regulated when one avoids double counting of glu- 
ons, thanks to a longitudinal cut-off. For example, let us take a cut-off fc+ in k + , 
and consider gluons with k + > fc+ only as part of the projectile, in order to regulate 
the integration with the weight (|10l) . and gluons with k + < fc+ only as part of the 
target semi-classical field, so that (5 XoXl ) and (S XqX2 S X2Xl } become fc+-dependent. 

The evolution with fc+ of \S X0X1 ), + , (5 XoX2 S X2X1 ),+ and similar expectation 
values is obtained by integrating out QCD quantum corrections slice by slice in k + . 
At leading log accuracy, one gets 

k c9 k +(S XaXl ) k + = I d 2 X 2 2 (5x0X2 5 X2 xi) fc + - (5 XoXl ) fe + (11) 

where xfj = (x,— x^) 2 , and similar other coupled equations forming Balitsky's 
infinite hierarchy^ involving expectations valu es of operators with an arbitrary 
number of Wilson lines. The JIMWLK equation^ 10 ! 11 ! 12 ! 13 ! 1 ^ ! is equivalent 
to this full hierarchy, and gives the evolution with fc+ of the target gluon field 
probability distribution (...) fc+ . The fact that the quantum corrections to the 
statistical distribution of gluons fields can be absorbed by evolution of the statistical 
distribution is a crucial self-consistency check for the CGC formalism. 

Alternatively, the same equation (fTT|) can be derived-^ by extracting the soft- 
gluon logarithmic divergence coming from the photon light-cone wave function in 
(j9|), and requiring the DIS cross section to be independent of the cut-off playing 
the role of a longitudinal factorization scale. The consistency of the two approaches 
shows the validity of the high-energy factorization formulae (j6]) and ([9|) for DIS 
once the cut-off fc+ is applied, since the B- JIMWLK evolution properly resums the 
high-energy leading logs. 

In the large N c limit, one has (5 Xt)X2 5 X2Xl ) fe + -4 (5 XoX2 ) fe + (5 X2Xl ) fc +, so that 

(fTTj) becomes a closed equation, called Balitsky-Kovchegov (BK) equation F^^^^ 
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When using the NLO factorization formula ©, by consistency, one should prob- 
ably resum not only the high-energy leading logs, but also the next-to-leading ones. 
This is performed by the NLL extension of the equation (ITT|) . derived recentlyJ^HH] 

Let us come back to the main assumption of the CGC formalism, i.e. that the 
target can be modeled by a random classical gluon field. First, we have seen the 
self-consistency of this assumption thanks to the JIMWLK equation. Second, when 
the gluon occupation number is large, it is large compared to commutators, so that 
the physics becomes effectively classical and the CGC assumption is justified. Due 
to Lorentz contraction, the valence partons of a relativistic large nucleus are densely 
packed, so that the CGC formalism applies to large nuclei as soon as relativistic 
effects are important. Moreover, the B- JIMWLK evolution leads to a more and 
more packed target, so that all this formalism is always valid at high enough energy, 
whatever is the target, even if in intermediate stages of the high-energy evolution, 
the target is dilute. Finally, the dilute limit (or two-gluon approximation) of the 
equation ([TTT) and of its NLL extension correctly reproduces the LL and NLL BFKL 
equations respectively, so that at this accuracy, no information about the dilute 
regime is lost when taking the semi-classical target approximation. However, this 
approximation for the target field might prevent to reproduce the full NNLL BFKL 
equation in the dilute limit. 

3. Large NLL corrections and collinear resummations 

The BFKL equation receives large corrections at NLLpHHH giving a pathological 
behavior to its solutions and signaling a breakdown of the corresponding formalism. 
Soon after this discovery, it has been realized that these issues come entirely from the 
collinear and anti-collinear kinematical regimes and not from the Regge kinematical 
regime in wh ich the B FKL equation is derived. The solution is then to perform 
resummationJ^^^HIH] j n collinear and anti-collinear regimes, in order to extend 
the validity of the NLL BFKL toward these regimes, which are in principle govern 
by DGLAP physics. Several successful resummations schemes have been proposed, 
which differ only by inessential higher order contributions. The general idea behind 
these schemes is to add appropriate NNLL and higher order terms to the NLL 
BFKL equation in such a way that the (N)LO DGLAP equation is reproduced in 
the collinear and anti-collinear limits. 

The physical origin of the large corrections has also been understood. They can 
be classified as follows. 

• Running coupling scale: The NLL BFKL kernel contains terms which are 
the beginning of the expansion of the running coupling. When the coupling 
is taken to run, but not with the correct scale, some of these terms stay in 
the NLL kernel instead of being resummed into the running coupling. If the 
running coupling scale differs in the collinear or anti-collinear regime from 
the one dictated by DGLAP physics, the left-over terms in the kernel are 
particularly large. 
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• Non-singular part of DGLAP splitting function: The LL BFKL equation 
contains the singular part of the LO DGLAP gluon splitting function both 
in the collinear and anti-collinear regimes, but not the non-singular part. 
The non-singular part shows up in an expanded way, term after term when 
going to NLL BFKL and further. This expanded form appears to be patho- 
logical, and one should instead include the full non-singular part together 
with the singular part of the splitting function in a modified version of the 
LL BFKL equation. 

• Kinematical constraint: The last large corrections, which are parametrically 
the largest ones, come from a bad treatment of kinematics when taking the 
high-energy limit too strictly. They can be resummed by imposing some 



collinear and/or anticollinear poles of the kernel in Mellin space!^ 

The corrections of the third type depend on the factorization scheme, in par- 
ticular on the precise definition of the evolution variable for the BFKL equation. 
They are however identical in any Yang-Mills theory. By contrast, the corrections of 
the first two types are high-energy factorization scheme independent. But they vary 
from theory to theory. In particular, they are both absent in Af = 4 SYM theory. 

Since the linearization of the NLL BK equation gives the NLL BFKL equation, 
those large corrections are also present in the NLL BK equation. The nonlinear- 
ity related to gluon saturation in the BK equation cannot help stabilize the solu- 
tions in th e pr esence of those large corrections, as shown by toy-model numerical 
simulations. Hence, it is necessary to perform collinear resummations also for the 
NLL BK equation before using it for phenomenology. However, the resummation 
schemes developed for NLL BFKL cannot be used directly for the NLL BK equa- 
tion, because of the nonlinearity of the latter and of the switch from transverse 
momentum to transverse position space. 

The terms related to running coupling in the NLL BK equation have been an- 
alyzed in detail, and an appropriate running coupling prescription is available 1381 
Concerning the non-singular part of the gluon splitting function, only a very rough 
treatment is available}^ so that more work on this issue is desirable. Let us now 
discuss the kinematical issues in the BK context. 

4. Resummation via kinematical improvement 

The need for the kinematical constraint in the evolution equation was first dis- 
covered in a covariant formalism, in a quite technical calculation aiming at the 
description of other effects like angular ordering!^ However, the origin of the kine- 
matical constraint and of the related large NLL corrections can be understood very 
clearly in light- front perturbation theory, as shown in Ref. 40 

The easiest derivation of the LL BFKL and BK equations is performed by calcu- 
lating the leading logarithmic contributions to the squared wave functions | ^ / qq g L . ( ^ ; | 
for the quark antiquark and n gluons Fock component of the projectile in the multi- 




BFKL equation, or by s hifting 
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Regge kinematics SHIHl In light-front perturbation theory, the kinematical issues are 
only due to the usual implementation of the multi-Regge kinematics on the energy 
denominators, coming with each additional gluon. In this context, the multi-Regge 
kinematics is the assumption that the k + of the partons are strongly ordered^ 

q+ > k+, q+-k+ » k+ » • ■ • » k+ +1 , (12) 

and that simultaneously their transverse momentum are of the same order 

Q 2 ~ k 2 ^ kl ~ • • • ~ k 2 +1 . (13) 

If one assumes only the k + ordering (IT21 . then for each diagram, the energy de- 
nominator associated with the emission of the gluon j, with 2 < j < n + 1, reduces 
to 

Q 2 q+k 2 k\ k 2 

— 1 - 1 — H — (14) 

2q+ 2k+(q+-k+) 2k+ 2k+' V ' 

Then, in the strict multi-Regge kinematics (|T2"j) - (fT3")) . one can approximate each 
energy denominator by its last term 

Q 2 g+ k 2 k| k| kj 

W + 2k+{q+-k+) + 2k+ + "' + 2kf-2kf~ 3 ■ (15) 

That approximation is crucial to the derivation of the LL BFKL and BK equations. 
When the high-energy factorization is formulated with a cut-off in k + , as presented 
in section[2l the fc + -ordering (| 12[) of the gluon cascades generated by the BFKL and 
BK equations is automatic. However, the integration over transverse momentum or 
position in the kernel of the LL BFKL and BK equations is not restricted, so that 
it contains contributions badly violating (|13[) and thus outside the multi-Regge 
kinematics. More precisely, once the ordering (|12l) is satisfied, the approximation 
(|T5j) is good typically for gluon cascades in the multi-Regge kinematics (fT3)l or in the 
anti-collinear kinematics, i.e. with k 2 growing with j, but not in the deep collinear 
regime, i.e. with k 2 decreasing fast enough with j. 

From this discussion, ones predicts that the fc + -ordered BFKL and BK equations 
are less reliable in the deep collinear regime, which is especially relevant for DIS, 
than in the rest of the available phase-space. This observation explains why the 
Mellin transform of the NLL kernel of the fc + -ordered BK equatiorPH contains a 
large and unphysical collinear triple pole, but no anti-collinear triple pole. 

All this discussion is however factorization scheme dependent. There are indeed 
other ways to regularize the rapidity divergence of Wilson line operators and the soft 
divergence (|10j) . For example, one can replace the cut-off in k + by a cut-off in fc~ . In 



a Here, the gluons are ordered along the gluon cascade from the closest to the projectile, 2, to the 
closest to the target, n + 1. is the momentum of the antiquark, and q+ the one of the virtual 
photon. 
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that case, the BFKL and BK equations implement automatically the /c~-ordering 

|^ + Wf^ <<fc2 " <0 "" <<fc " +1 ' (16) 
so that approximations of the type (|15l) are good. However, this time, the approx- 
imation of the energy denominators by the expression (fT4|) is not always valid. In 
fact, we are in the inverse situation: the approximation of energy denominators by 
the k~ of the last gluon is in general valid except in the deep anti-collinear regime. 
Hence, a spurious anti-collinear triple pole is expected in the Mellin transform of 
the fc~-ordered NLL BFKL and BK kernels, but no collinear triple pole. 

In fact, whatever factorization scheme is chosen, the approximation of the energy 
denominators by the k~ of the last gluon is valid if and only if both the k + and k~- 
ordering flT2"j) - flTBl) are satisfied, whereas the assumption (ITB")) is unnecessary. More- 
over, one can show that if either (fl2l) or (|16l) is not satisfied, no leading logarithmic 
contribution is generated!^ Hence, the BFKL and BK equatio ns sho uld actually 
resum precisely the gluon cascades ordered both in k + and fc~ | 40 | 42 ] wnereas the 
usual assumption (|1 3|) is unnecessary and misleading. In any factorization scheme, 
imposing both k + and k~ ordering allows to resum into a modified LL equation the 
largest unphysical collinear or anti-collinear contributions appearing in the kernel 
at higher logarithmic orders, corresponding in Mellin space to triple poles at NLL, 
quintuple poles at NNLL and so on. 

Only the fc + -ordcring or the fc~-ordering, not both, can be made automatic 
by a choice of high-energy factorization scheme. The other one has to be included 
as a modification of the LL kernel, by restricting the integration over transverse 
momentum or position to the domain allowing to satisfy both (p~2]) and (fTH]). For 
the BFKL equation in momentum space, this restriction of the integration domain 
is precisely the kincmatical constraint found in Ref. 1331 

So far, we have shown the necessity of the kinematical constraint in the BFKL 
and BK evolutions by an analysis in transverse momentum space. However, evo- 
lutions equations including gluon saturation like the BK and JIMWLK equations 
are more easily written in transverse position space. Hence, one has to translate 
the kinematical constraint to a restriction in position space. A first implementation 
of the kinematical constraint in position space for the BFKL and BK equations 
has been proposed in Ref. |40] I have recently clarified^ some issues related to this 
first implementation, in particular by doing a more careful treatment of virtual 
corrections and by analyzing systematically the factorization scheme dependence of 
the kinematical constraint. In the end, one gets a modified BK evolution equation 
with memory effects and a restricted phase-spa ce f or gluon emission which gets 
progressively wider in the course of the evolution!^ 
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